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CHAPTER 1 


INTRODUCTION 

This report describes work carried out by the Flight Research 
Laboratory of the University of Kansas (KU-FRL) and sponsored by 
Grant NSG 1574 from the National Aeronautics and Space Administration 
(NASA), Langley Research Center. The purpose of this project was 
to correlate theoretically predicted aerodynamic characteristics 
of the Ayres Thrush S2R-800 airplane with full scale wind tunnel 
data. The theoretical prediction schemes have been documented in 
this report to provide an analytical method for calculating selected 
longitudinal and lateral-directional characteristics of the Ayres 
Thrush, winglets on and winglets off. 

It has been previously found that winglets have favorable 
effects on wing tip vortex entrainment and wake interaction charac- 
teristics of chemical sprays dispensed from the airplane. These 
effects are especially Important for an agricultural type airplane 
such as the S2R-800; therefore, this report investigates the wing- 
let effect on the S2R-800 lateral-directional stability. 

The two methods used in this report for theoretical predictions 
were a Quasi-Vortex Lattice Method (QVLM, Reference 3), which is a 
computer method, and DATCOM analytical methods (Reference 1). The 
results of these calculations are compared to full scale wind tunnel 
data where possible, and recommendations and conclusions given con- 
cerning these comparisons. All calculations have been done in the 
stability axes system. 
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CHAPTER 2 


THE AYRES THRUSH S2R-800 AIRPLANE 

The Ayres Thrush S2R-800 is designed especially for agricultural 
flying. Typical applications for the S2R-800 include seeding, 
fertilization, insect control, and defoliation. The S2R-800 is a 
monoplane featuring a full cantilever low wing and, with the exception 
of fabric empennage surfaces, is of all-metal construction. The 
major physical characteristics are listed in Table 2.1, and a three- 
view drawing is shown in Figure 2.1 

The winglets used on the S2R-800 consist of a modified GA(W)-2 
airfoil. Table 2.2 lists the airfoil coordinates of the winglets. 
Figure 2.2 shows a three-view sketch of the aircraft with winglet 
installation. 
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Table 2.1: Specifications of the Ayres Thrush S2R-800 


Wing : 

Airfoil 

Section Characteristics 
deg 



, per deg 
a 

deg 

Area, (ft^) 

Exposed area, (ft^) 

Span, m (ft) 

Exposed Span, m (ft) 

Aspect ratio 
Exposed aspect ratio 
Thickness ratio 
Dihedral , deg 
Taper ratio 
Root chord, m (ft) 

Mean aerodynamic chord, m (ft) 

Incidence angle at root, deg 

Incidence angle at tip, deg 

Sweep angle of quarter chord line, deg 

Horizontal tail : 

Airfoil 

Area (including elevator and tabs), ra^ (ft^) 

* 

See Section 2.1, page 2.8 and Figure 2.10. 


NACA 4412 

0.4933 

0.5067 

0.105 

7.5 

30.34 (326.6) 
27.45 (295.5) 
13.27 (43.5)* 
12.00 (39.4) 
5.81 
5.25 
0.12 

3.5 
1.0 

2.29 (7.5) 
2.29 (7.5) 

0 

-1.5 

0 

NACA 0003 
5.25 (56.47) 
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Table 2.1; (continued) 

Span, m (ft) 

Aspect ratio 
Thickness ratio 
Dihedral , deg 
Taper ratio 
Root chord, m (ft) 

Mean aerodynamic chord, m (ft) 

Incidence angle at root and tip, deg 
Quarter chord sweep angle, deg 

Vertical Tail ; 

Airfoil 

Area (including rudder) , (ft^) 

Span, m (ft) 

Aspect ratio 
Taper ratio 
Root chord, m (ft) 

Mean aerodjmamic chord, m (ft) 

Quarter chord sweep angle, deg 

Winglet ; 

Airfoil 

Area, (ft^) (per winglet) 

Span, m (ft) 

Aspect ratio 
Taper ratio 


4.71 (15.45) 

4.23 

0.03 

0 

0.66 

1.34 (4.4) 

1.13 (3.71) 

0 

8.5 

NACA 0003 
2.12 (22.77) 

1.51 (4.96) 

1.08 

0.5 

1.87 (6,14) 

1.3 (4.26) 

14.0 

Modified GA(W)-2 
1.80 (19.34) 

1.52 (4.98) 

1.28 

0.56 
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Table 2.1: (concluded) 


Root chord, m (ft) 1.52 (4.98) 

Mean aerodynamic chord, m (ft) 1.22 (3.99) 

Leading edge sweep angle, deg 12.5 

Weights : 

Typical operating weight, N (lb) 

Empty weight, N (lb) 

Powerplant : 

Number of engines 
Manufacturer 
Model 

Takeoff power, kw (hp) 

Takeoff rpm 


Hamilton Standard 
3 D40/EAC 
3 

2.74 (9) 


Propeller : 

Manufacturer 

Model 

Number of blades 
Diameter, m (ft) 


1 

Wright 

R-1300-1B Cyclone 
1072.4 (800) 

2600 


34,696 (7800) 
18,238 (4100) 


2.4 



Table 2.2: Airfoil Coordinates for Winglets 


z/c for - 



Upper surface 

0 

0 

.0020 

.0077 

.0050 

.0119 

.0125 

.0179 

.0250 

.0249 

.0375 

.0296 

.0500 

.0333 

.0750 

.0389 

.1000 

.0433 

.1250 

.0469 

.1500 

.0499 

.1750 

.0525 

.2000 

.0547 

.2500 

.0581 

.3000 

.0605 

.3500 

.0621 

.4000 

.0628 

.4500 

.0627 

.5000 

.0618 

.5500 

.0599 

.5750 

.0587 

.6000 

-.0572 

.6250 

.0554 

.6500 

.0533 

.6750 ' 

.0508 

.7000 

.0481 

.7250 

.0451 

.7500 

.0419 

.7750 

.0384 

.8000 

.0349 

.8250 

.0311 

.8500 

.0270 

.8750 

.0228 

.9000 

.0184 

.9250 

.0138 

.9500 

.0089 

.9750 

.0038 

1.0000 

-.0020 


Lower surface 
0 

-.0032 
-.0041 
-.0060 
-.0077 
-. 0090 . 
-.0100 
-.0118 
-.0132 
-.0144 
-.0154 
-.0161 
-.0167 
-.0175 
-.0176 
• -.0174 
-.0168 
-.0158 
-.0144 
-.0122 
-.0106 
-.0090 
-.0071 
-.0052 
-.0033 
-.0015 
.0004 
.0020 
.0036 
.0049 
.0060 
.0065 
.0064 
.0059 
.0045 
.0021 
-.0013 
-.0067 
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Figure 2.2: Three-view sketch of Ayres Thrush S2R-800 showing winglet installation (Uefercuice 5) 



2.1 Geometric Parameters of Fuselage. Wing, Wlnglet. and Tails 

Many of the geometric parameters used in the subsequent analysis 
were taken from manufacturer’s specifications or actual blueprints. 
Pertinent dimensions for the fuselage, wing, wlnglet, horizontal 
tail, and vertical tall are shown in Figures 2.1.1, 2.1.2, 2.1.3, 
and 2.1.4, respectively. It should be noticed that the wing span as 
defined in Figure 2.1.2 is different from the wing span definition 
in Figure 2,1, This is done for convenience in the analysis so that 
the wing may be considered to be an untapered wing, i,e,, the tapered 
wing tips have been replaced with an untapered tip of equivalent area. 
This has little, if any, significant effects on the accuracy of the 
analysis as evidenced by Figures 3,1.1 - 3,2,3* Care must be taken 
not to generalize this conclusion to other airplane configurations. 
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Figure 2.1.2: Definition sketch of wing dimensions 
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2.2 Flight Condition Used in the Analysis 

The center of gravity of the airplane was fixed at 25 percent 
of the wing mean aerodynamic chord in the longitudinal direction and 
at 11 percent of the wing mac above the wing root. See Figure 2.1. 

Since the S2R-800 is a low subsonic regime vehicle, the flow has 
been assumed to be incompressible; therefore, the results of the analysis 
are valid for both the climb and cruise speeds of the Ayres Thrush. 
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CHAPTER 3 


PRESENTATION OF RESULTS 

In this chapter the results of the analysis are presented and 
discussed. Full scale wind-tunnel data (q = 718.2 N/m^, Reference 4) 
are compared to theoretical results where applicable. These results 
apply to both cruise and climb speeds of the Ayres Thrush S2R-800. 

3.1 Lift Curve 

The analytical method used in predicting the lift curve of the 
S2R-800 is presented in Chapter 4. The Quasi-Vortex-Lattice Method 
(QVLM) prediction of the lift curve is discussed in Section 6.2. 

The results of both methods are compared to each other and to full- 
scale wind-tunnel data in Figure 3.1.1. 

Figure 3.1.1 shows that both predictions agree very well with 
wind tunnel data. The QVLM prediction underestimates the lift curve 
slope because QVLM assumes that the wing dominates the lift behavior 
of the total airplane and therefore doesn’t take into account fuselage 
and/or empennage lift effects. 

3.2 Sideslip Derivatives 

In Figures 3.2.1 through 3.2.3, the calculated sideslip deriva- 
tives for the winglet off configuration are compared to full scale 
wind tunnel data. The calculations compare favorably with the tunnel 

data. In Figure 3.2.2 the predicted C increases with increasing 

8 

angle of attack, while the tunnel data shows C decreasing with 

6 

increasing angle of attack. This is to be expected, since the 
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analytical method taken from Reference 1 does not account for 


fuselage and wing-fuselage Interference effects on C 


n 


However, 


6 


the nominal predicted value of C agrees well with the average 

”6 

tunnel C . 

6 

Comparing Figures 3.2.4 - 3.2.6 to 3.2.7 - 3.2.9, it can be 
seen that wlnglet effect on the sideslip derivatives has not been 
properly accounted for in the analytical method of Reference 1. 

This is especially evident in the sideslip derivative C. , where 
the computer results (Reference 3) indicate that winglets have a 
very strong influence on C. , while the analytical results (Reference 1 ) 
show only a weak influence. 

Figures 3.2.7 - 3.2.9 show that winglet cant angle has a 
significant influence on the sideslip derivatives. No wind tunnel 
data for the airplane-winglet configuration were available. 
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Figtire 3.1.1: Lift 
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Figure 3«2.1: Comparison of predicted C , winglets off, 
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based on QVLM computation for the 
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CHAPTER 4 


PREDICTION OF PROPELLER-OFF, HORIZONTAL TAIL-OFF LIFT CHARACTERISTICS 
~ USING AN ANALYTICAL METHOD 

In this chapter the propeller-off, horizontal tail-off lift behav- 
ior of the Ayres Thrush will be discussed. This chapter presents the 
methods used in calculating the lift characteristics of the Thrush and 
illustrates the contribution of each relevant component to these char- 
acteristics. Table 2.1 lists the pertinent parameters and their magni- 
tudes used in the analysis. All calculations were made in English 
units. Extensive use of Reference 1 was made to predict the lift be- 
havior of the Thrush. 

4.1 Zero Lift Angle of Attack, 

The zero lift angle of attack at low speeds of the complete air- 
plane, minus the horizontal tail, is considered to be relatively inde- 
pendent of fuselage effects and is primarily determined by the wing 
airfoil properties. Therefore, the airplane zero lift angle will be 
found by considering only the wing including the effect of wing twist. 

4.1.1 Wing Contribution, 

For untwisted constant section wings the zero lift angle of 
attack is given by Equation 4. 1.1.1 from Reference 1. 

(a)g^Q “ “i “ c“ (4. 1.1.1) 

a 

where: 

is the angle of attack for the xd.ng section design lift 
coefficient, obtained from Table 2.1 
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C, is the wing section design lift coefficient, from Table 2.i 

Is the wing section lift curve slope, from Table 2.1 
a 

To account for wing twist. Equation 4. 1.1. 2 from Reference 1 
Is used. 

Aa 

® (4. 1.1.2) 


where . 

Ao 

— g— Is the change In wing zero-lift angle of attack due to a 
unit change In linear wing twist, obtained from Figure 4. 1.1.1 

0 Is the twist of the wing tip with respect to the root section. 
In degrees (negative for washout), from Table 2.1. 

Table 4. 1.1.1 summarizes the wing contribution to the airplane 
zero lift angle of attack. 


4.1.2 a of the Avres Thrush 

o 

The zero lift angle of attack of the complete airplane, prop-off, 
horizontal tall-off. Is considered to be Identical to the wing-alone 
a^. Therefore, for the Ayres Thrush: 

= -3.69 deg 

4.2 Llft-Cuirve Slope, 

a 

For a conventional horizontal tail-off configuration, the lift- 
curve slope in the linear angle-of-attack range can be found by con- 
sidering the following components: 

(1) Wing, including Interference effects 

(2) Fuselage, including nose lift and fuselage-wing interference. 
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4.2.1 Wing Contribution to 

ct 

The wing-alone lift-curve slope, based on exposed planform 
dimensions, for a straight tapered wing in the low subsonic region 
can be calculated using the standard Polhamus equation. 


a (e) 


2TTAg/57.3 


■Ji 


per deg 


2 +,/ — (1 + tan^ A ) + 4 


w 


c/2 


(4. 2. 1.1) 


where: 

A is the aspect ratio of the wing based on its exposed planform, 
e 

obtained from Table 2.1. 

C„ 


K = 


a 


2ir 


A / is the wing sweep angle at the half-chord location, from 
c/2 

Table 2.1. 

Equation 4. 2. 1.2 from Reference 1 accounts for wing-fuselage 
interference on the wing contribution to . 


(C. ). 


(C, ), 


-a "(f) "(f) "(e) 


per deg 


(4. 2. 1.2) 


w 


where! 

K is the ratio of the wing lift in the presence of the fuse- 

^(f) 

lage to the wing-alone lift, obtained from Figure 4. 2. 1.1. 

S is the exposed wing area, obtained from Table 2.1. 
e 

The summary calculation for the wing contribution to is 

a 

shown in Table 4*2.1«1* 

For the Ayres Thrush: 

(C, ) = 0.0732 per deg 

a (f ) 
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4.2.2 Fuselage Contribution to Cj^ 

a 

The fuselage contribution to is accounted for in Equation 

a 

4. 2. 2.1 taken from Reference 1. 


a (w) 


>w, , f 

(w) a (e) W 


per deg (4. 2. 2.1) 


where • 

K-. is the ratio of the fuselage lift in the presence of the 
(w) 

wing to the wing alone lift, obtained from Figure 4. 2. 1.1. 

tr 2 

Kjj = stTJTc — ) ^"s — ^ fuselage nose lift based on 

a (f) 

slender body theory where: 

r is the radius of the equivalent circular fuselage, obtained 
from Figure 2.1.1. 

The summary calculation for the fuselage contribution to 

is shown in Table 4. 2. 2.1 

For the Ayres Thrush: 

(C^ ') f. = 0.0130 per deg. 
a (w) 

4.2*3 C- of the Complete Airplane, Horizontal Tall Off 


a 


The horizontal tail-off, power-off lift-curve slope of the 
airplane is given as: 


^w_, h, , P®'' (4. 2. 3.1) 

a a (f) a (w) 

Table 4. 2. 3.1 summarizes the contribution of each component. 

For the Ayres Thrush: 

Cj^ = 0.0862 per deg. 

a 
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4.3 Lift-Curve in the Non-Linear Region 


To obtain a description of the non-linear portion of the lift 
curve, propeller and horizontal tail off, it will be necessary to 
consider the following components: 

(1) the wing and wing-fuselage 

max 

(2) the wing and wing-fuselage angle of attack at 

max 

(3) ■ the upper linearity of the wing lift-curve slope. 


4.3.1 of the Wing, Wing-Fuselage 

max 

The TnaviTnntn lift Coefficient of a wing with twist may be estimated 

from the assumption that is reached when the local section lift 

max 

coefficient, C., at any position along the span is equal to the local 

A* 

C for the corresponding section. The method used is taken from 
max 

Reference 2. 

The first step in finding of the wing is to calculate the 

max 

variation of the local section lift-coefficient with span location, 
at a total C. of 1. This is done with Equation 4. 3. 1.1, which only 

Li 

applies to unswept, untapered, linearly twisted wings. 

(4. 3. 1.1) 

a b 

where ; 

C is the wing section lift coefficient due to wing angle of 
a 

attack, given by Equation 4. 3. 1.2 

C is the wing section lift coefficient due to wing twist, given 
^b 

by Equation 4. 3. 1.3. 
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(4. 3. 1.2) 


a 


Aa 

+ -f> 


(4. 3. 1.3) 


a a 

where; 

^1* ^2* ^3’ ^4 additional and basic lift 

distributions, obtained from Figure 4. 3. 1.1 

n = and is the wing spanwise station 

9 is the wing twist measured from the wing root in degrees, 
negative for washout, from Table 2.1 
Aa 

is the ratio of the change in wing zero lift angle of attack 
with wing twist, obtained from Figure 4. 1.1.1. 

Table 4. 3. 1.1 summarizes the calculation of the wing lift distri- 
bution. From this table the minimum value of the ratio of (C. - C, ) 

x» x», 

max b 

to @ ~ 1 is considered to be the maximum lift coefficient of the 

a 

wing; is the wing section maximum lift coefficient. Table 4. 3. 1.1 

max 

also summarizes the calculation of Cj^ for the wing. 

max 

For the Ayres Thrush: 

)„ » 1.412 

max 

To account for the presence of the fuselage. Reference 1. gives the 
following equation: 


^wf 

max 


«=L >wf 
max . 

(C, ) 'w 

L w max 


). 


(4. 3. 1.4) 


max 


where! 


("l >wf 

max 

max 


is the ratio of the wing and fuselage C 


to the wing- 


max 
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alone , obtained from Figure 4. 3. 1.2. 

max 

Table 4. 3.1.2 summarizes the calculation for of the wing- 

max 

fuselage. 

For the Ayres Thrush: 

(C^ ),f= 1-412. 

max 


4.3.2 The Wing and Wing-Fuselage Angle of Attack at 

max 

For high-aspect-ratio, constant-section wings, the angle of attack 
at the wing is computed using Equation 4. 2. 3.1 from Reference 1. 


(Cl >w 

Cl ^ (Cl ° ^L 

max a (e) “*ax 


(4. 3. 2.1) 


where : 


(C, ) is obtained from Section 4.3.1 and is the wing maximum 

L w 
max 

lift coefficient 

(C, ) is obtained from Section 4.2.1 and is the exposed wing 
^o.-'Ce) 

lift-curve slope 

a is obtained from Section 4.1.1 and is the wing zero lift 
o 

angle of attack 

Aa is the angle of attack increment for subsonic maximum 

max 

lift, obtained from Figure 4. 3. 2.1. 

To account for the fuselage. Reference 1. gives the following 
equation : 


(“Cj^ ^wf 

^ X max / N 

max L max 

max 


(4. 3. 2. 2) 


4.7 



where ; 


r — Is the ratio of wing-fuselage angle of attack at 

w max 

max 


to the wing-alone angle of attack at , obtained from Figure 4. 3. 2. 2, 


max 


A summary calculation for the angle of attack at is shown 

Li 


max 


in Table 4.3.2 .1. 

For the Ayres Thrush: 

*C. 


(“r >wf **17.2 deg. 


max 


4.3.3 Upper Limit of Linearity of the Lift Curve 


The angle at which the lift-curve slope is no longer linear, 
for the tail-off Thrush configuration, is considered to be approxi- 
mately equal to the corresponding angle for the wing-alone configu- 
ration. From Reference 1: 

+ + 

a * + ~f~ • 0 deg (4. 3. 3.1) 

where : 

is the section angle of attack, in degrees, at which the lift 
curve slope is no longer linear, from Table 2,1. 

A summary calculation is shown in Table 4. 3. 3.1. 

For the Ayres Thrush: 

o"** => 8.1 degrees 

4.4 Summary 


In this chapter an analytical method for predicting the lift behav- 
ior of the Ayres Thrush was presented. Table 4.4.1 below lists the 
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pertinent lift characteristics for the Ayres Thrush, horizontal tail 
and propeller off. Figure 4.4.1 compares these predictions to actual 
full-scale wind tunnel data. 


1 Table 

4.4.1: Lift Characteristics of 

Tail and Propeller Off 

the Ayres Thrush 

» 

Symbol 

Description 

Reference 

Magnitude 

a 

o 

Zero lift angle of attack, 
deg. 

Section 4.1.1 

-3.69 

a 

Linear lift-curve slope, 
per deg. 

Section 4.2.1 

0.0862 

max 

Maximum lift coefficient 

Section 4.3.1 

1.412 

1 \ 

Angle of attack at the 
maximum lift coefficient. 

Section 4.3.2 

17.2 

» max 

deg. 



1 

Angle of attack for lift- 

Section 4.3.3 

8.1 


curve slope is no longer 
linear, deg. 
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Table 4. 1.1.1: Wing Contribution 

to a 

0 


Symbol 

Description 

Reference 

Magnitude 

“i 

Angle of attack at wing section design 
lift coefficient, deg. 

Table 2.1 

0.4933 


Wing section design lift coefficient 

Table 2.1 

0.5067 

a 

Wing section angle of attack, per deg. 

Table 2.1 

0.105 

p 

0 

CD 

II 

0 

Zero lift angle of attack for untwisted 
wing, deg. 

Equation 4.1il.l 

-4.33 

Aa 

r\ 




U 

0 

Change in zero lift angle of attack due 
to wing twist 

Figure 4. 1.1.1 

-0.427 

0 

Wing twist, negative for washout, deg. 

Table 2.1 

-1.5 

<“o >0 

Zero lift angle of attack for twisted 
wing, deg. 

Equation 4. 1.1. 2 

0.6405 


Summary: («q) 0 ^q + (<* 0)0 “ -3.69 deg. 


i 


Symbol 


w 


c/2 


a (e) 


w 


(f) 


w 


Table 4 . 2 . 1 . 1 : Wing Contribution 

to 

a 


Description 

Reference 

Magnitude 

Exposed wing aspect ratio 

Table 2.1 

5.253 

C /2tt 


0.9576 

Wing sweep at half chord, deg. 

Table 2.1 

0 

Exposed wing lift-curve slope, per deg. 

Equation 4. 2. 1.1 

0.0735 

Wing-fuselage interference factor 

Figure 4. 2. 1.1 

1.1 

Exposed wing area, (ft^) 

Table 2.1 

27.45 (295.5) 

Total wing area, (ft^) 

Table 2.1 

30.34 (326.6) 


(Ct )„ = 0.0732 per deg. 

a ”(f) 


Summary: 



Table 4.2. 2.1: Fuselage Contribution to C- 

L 

a 


Symbol 

Description 

Reference 

Magnitude 


Fuselage-wing interference factor 

Figure 4. 2.3.1 

0.16 

(w) 




Nose lift factor 

Section 4.2.2 

0.0359 

a (e) 

Wing lift-curve slope based on the 
exposed wing geometry, per deg. 

Section 4.2.1 

0.0735 

S 

e 

Exposed wing area, (ft^) 

Table 2.1 

27.45 (295.5) 

S 

w 

Total wing area, (ft^) 

Table 2.1 

30.34 (326.6) 


Siinunary: (C. = 0.0130 per deg. 

a ^w) 



Table 4. 2. 3.1: Linear of the Ayres Thrush 

a 







^prop off 
° horiz. tail off 


Table 4. 2. 1.1 


Table 4. 2. 2.1 


Equation 4 .2. 3.1 


.0732 


0.0130 


.0862 


Table 4. 3. 1.1: Maximum Wing Lift Coefficient 



C 


il 

a 



Eg. 4.3.1. 2 Eg. 4. 3. 1.3 



0 . 

1.160 

.0322 

1.412 

0.1 

1.156 

.0246 

1.424 

0.2 

1.543 

.0170 

1.447 

0.3 

1.120 

.00941 

1.482 

0.35 

1.106 

.00572 

1.505 

0.4 

1.088 

.00211 

1.533 

0.45 

1.168 

- .00138 

1.565 

0.5 

1.045 

- .00472 

1.603 

0.6 

0.988 

- .0108 

1.701 

0.7 

0.914 

- .0159 

1.844 

0.8 

0.816 

- .0195 

2.070 

0.9 

0.675 

- .0205 

2.505 

1.0 

0.300 

- .0110 

5.603 


4.13 



4.14 


Table 4. 3. 1.2: Summary Calculation for Wing-Fuselage C. 


max 


Symbol 


Description 


Reference 


Magnitude 


^wf 

max 

max 


Ratio of the wing and fuselage to 

max 

the wing-alone 

max 


Figure 4. 3. 1.2 1.0 


Summary: (Cj^ = 1.412 

max 



4.15 


• > 





Table 4. 3. 2.1: Summary Calculation 

for a 

max 


Symbol 

Description 

Reference 

Magnitude 

«L )„ 

max 

Wing maximum lift coefficient 

Section 4.3.1 

1.412 

a (e) 

Exposed wing lift-curve slope, per deg 

Section 4.2.1 

.0735 

a 

o 

w 

Wing zero lift angle of attack, deg. 

Section 4.1.1 

-3.69 

max 

Increment in angle of attack for wing 
maximum lift, deg. 

Figure 4. 3. 2.1 

1.2 

<“Cj^ >wE 
max 

(“Cl >w 

max 

Ratio of wing-fuselage angle of attack 
at Cj^ to the wing-alone angle of 

max 

attack at 

max 

Figure 4. 3. 2. 2 

1.03 

Summary: (a^ “ 17.2 deg. 

^max 





Table 4. 3. 3.1: 


Summary for a 


Calculations 


Symbol 

Description 

Reference 

Magnitude 

+ 

“l 

Section angle of attack at which the lift- 
curve slope Is no longer linear, deg. 

Table 2,1 

7.5 

Aa 

6 

Incremental zero lift angle due to wing 
twist 

Figure 4.1. 1.1 

-0.43 

0 

Wlng-tlp twist with respect to the wing 
root, neg for washout, deg. 

Table 2.1 

-1.5 


Summary : 


+ 

a 


8.1 deg. 



3 


QUARTER-CHORD SWEEP ANGLE, (deg) 

Figure 4. 1.1.1: Change in Wing Zero Lift Angle of Attack 

Due to Wing Twist (Reference 1) 






body diameter , ^ 
wing span b 


Figure 4. 2. 1.1: 


Lift Ratios K ,c\ and . 

w(f) f(w) 

(Reference 1) 


4.18 















1.3 



b 

Figure 4. 3. 2. 2: Angle of Attack for Maximum Lift. 

(Reference 1) 


.22 


prop off 
wit off 
tail off 

per deg 



a - deg 

Figure U.U.l: S2R-800 lift curve based on DATCOM methods. 
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CHAPTER 5 


PREDICTION OF PROPELLER-OFF LATERAL-DIRECTIONAL 
STATIC STABILITY USING AN AYALYTICAL METHOD 


In this chapter an analytical method for predicting the 
propeller-off, lateral-directional behavior of the Ayres Thrush 
will be discussed. The derivatives that are- considered here 


are the side force due to sideslip derivative, C ; dihedral effect, 

C ; and directional stability, C . The methods used are primarily 
^8 '^8 
from Reference 1. 


5.1 Side Force Derivative, C 
^8 

The side force due to sideslip, C , of the complete airplane 

^8 

is found by considering the contributions of the following components: 

(1) Wing 

(2) Fuselage 

(3) Vertical tail 

(4) Winglets 

Unless the horizontal tail has large twist or dihedral, it can be 
safely ignored in the calculations. These contributions to C can 
be represented by: 

^*^y ^PROP OFF “ ^So^v(wfh) C5.1.1) 

^ WLT OFF ® ^ ^ ® 


^Sg^PROP OFF “ ^^y ^PROP OFF ^^y ^WLT 
^ WLT ON ^ T^LT OFF ^ 


(5.1.2) 


5.1.1 Wing Contribution 


The wing contribution to C is primarily due to wing dihedral. 

^8 



This be computed by Equation from Reference 1. 

(C ) = -0.0001(r„) per deg (5. 1.1.1) 

>'6 “r 

where ; 

r is the wing dihedral angle in degrees, 
w 

For the Ayres Thrush the wing dihedral angle is 3.5 degrees. 

Therefore, 

(C ) = -.00035 per deg 

5.1.2 Fuselage Contribution to C 

The fuselage side-force due to sideslip contribution can be 
considered as the sum of the side forces on body and the wing- 
body interference. The fuselage alone is the main contributor. 

The wing-fuselage interference is primarily a function of wing 
vertical position on the fuselage. The total fuselage contribution 
to C at subsonic Mach numbers is given by Equation 5. 1.2.1 from 
Reference 1. 

_2/3 

(C ). = K. (C^ )_ V- per deg (5. 1.2.1) 

where : 

is the wing-fuselage interference factor, obtained from 
Figure 5. 1.2.1 

(C ) is the body alone side force due to sideslip. For an 
estimation of the body side force diie to sideslip, slender-body 
theory can be used, which gives (Cy. )g ” —0.0195 per deg. 

a 

V is the fuselage volume, obtained from Figure 2.1.1 
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Is Che wing area, from Table 2.1 
Table 5. 1.2.1 Is a summary calculation for the fuselage 
contribution to C 

^3 


5.1.3 Vertical Tall Contribution to C 


The vertical tail contribution to C is affected by the location 

^3 

of the horizontal tall, the fuselage crossflow on the vertical tail, 
and the wing-fuselage- induced sidewash. 

Reference 1 accounts for horizontal tail and crossflow effects 
by computing an effective aspect ratio. The vertical tail effective 
aspect ratio is: 



v^ A 


1 1 + K 

I h 


r ^v(fh) 

A 

. V(f) 



(5. 1.3.1) 


where? 


v(f) 

— — is the ratio of the aspect ratio of the vertical tail 

V 

in the presence of the fuselage to that of the isolated tail, 
obtained from Figure 5. 1.3.1 

is the geometric aspect ratio of the vertical tail, from 
Table 2.1 


ratio of the vertical tail aspect ratio in the 
presence of the horizontal tail and fuselage to that of the panel 
in the presence of the fuselage alone, obtained from Figure 5. 1.3.1 

iS- a factor accounting for the relative size of the horizontal 
and vertical tails, obtained from Figure 5. 1.3,1, 

Table 5. 1.3.1 shows the summary calculations made to obtain the 
effective aspect ratio of the vertical tail. 
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The effective aspect ratio found is used to calculate the 
lift-curve slope of the vertical tail. The standard Polhamus 
equation is used for this calculation. 


2ir A„ 


a ^(fh) 


eff 


per rad 



2+,/— T" (1 + tan2(A^/2^ )+ 4 


(5. 1.3.1) 


where: 


K 


a 


V 2tt 

(A /«) Is the mid-chord sweep angle of the vertical tail, 

^ c/2'v 

obtained from Table 2,1 

Table 5. 1.3.1 summarizes the calculation for computing (C^ )^. 

a 

The complete vertical tail contribution to C is given in 

. P ■ 

Equation 5.1. 3.2 from Reference 1. This equation adjusts the vertical 
tail lift-curve slope to account for wake and sidewash effects. 


■So^v(wfh) “ "4^\^v(fh) 3B^ - s per rad 

P « ^oo 


(5. 1.3. 2) 


where: 

is a factor which accounts for the relative size of the 
fuselage near the vertical tail to the size of the tail, from 

Figure 5^ 1 . 3 . 1 . 


30v St' 


S /S 
V w 


0.4 Z 


.725 3.06 


w 


+ 0.009 A 


cM^v ""f'w 


w 


(5. 1.3. 3) 


where: 


(A ,.) is the quarter-chord sweep of the vertical tail, obtained 
c/4 V 


from Table 2.1 
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is the vertical distance from the centerline of the equivalent 

fuselage to the quarter-chord point of the root chord of the exposed 

wing panel, obtained from Figure 2.1.1 

(w^)^ is the width of the equivalent circular fuselage at the 

wing, obtained from Figure 2.1.1. 

Table 5. 1.3.1 shows the summary calculations used to find the 

vertical tail contribution to C . 

^6 


5.1.4 Winglet Contribution to C 

Strictly speaking, at the present time there exist no analytical 
methods for calculating winglet effect on C . However, in this 
section this contribution will be approximated by treating the 
wlnglets as vertical tails at the wing tips. 

Twin vertical tails are treated in Reference 1. This method 
includes the effect of sidewash. Depending on winglet geometry, 
an effective aspect ratio is calculated from Figure 5. 1.4.1: 


^Sft^WLT 

^ r =0 

WLT 


"^S-^v(wfh) 

(Cy )v 

^3 eff • 


2 S 


(C ) 


WLT 


y-'v -- 57.3 S 
•^3 eff w 


per deg 


(5. 1.4.1) 


where : 

^^yg^v(wfh) 

- 7 -z — r is a mutual interference factor, obtained from Figure 

'■S ‘'v 
^3 eff 

5.1.4. 1 

(C ) is the lift-curve slope of one vertical-tail panel, 

^3 ^eff 

per rad, obtained from Figure 5. 1.4.1. 

To account for winglet cant angle effect on the side-force 
derivative, the following expression is used: 


5.5 



, per deg 


(5. 1.4. 2) 


2S 




WLT 


WLT 


-.0001 5 ^ 


where: 


r is the winglet cant angle 
WLT 

Summary calculations are shown 


(see Figure 2.2). 
in Table 5. 1.4.1. 


5.1.5 C of the Complete Airplane, Winglets On and Winglets Off 

The side-force due to sideslip derivative, , of the Ayres 

P 

Thrush, winglets and power off, is; 

(C ) “ -.0056 per deg. 

'■'^y ''prop off 

^ 1-ILT OFF 

Table 5. 1.5.1 summarizes the calculations and lists the effect of 
winglet cant angle on of the airplane. 



5.2 Yawing Moment Derivative, C 

8 

The weathercock stability, C , Is found by considering the 

6 

contributions of the following airplane components: 

(1) Wing 

(2) Fuselage, Including wing-fuselage Interference 

(3) Tails 

(4) Wlnglets 

These contributions to C can be represented by : 

“b 

^^n ^PROP OFF “ ^V^f(w) "^^^no ^v(wfh) (5.2.1) 

^ WLT OFF ^ ^ ^ 

^^n ^PROP OFF “ ^^n ^PROP OFF '** ^^n ^WLT (5.2.2) 

^ WLT ON ^ WLT OFF ^ 


5.2.1 Wing Contribution to C 

"b 


The wing contribution to weathercock stability Is primarily due 
to the asymmetrically Induced drag distribution caused by an asym- 
metrical lift distribution. 

For low subsonic speeds. Reference 1 gives the yawing moment 
derivative as; 

f , tan A ,, 

(C ) w 4 1 c/4 

n gW 57.3 I 4'irA irA (A +4 cos A ,,) 

B I w w' w c/4 


[ 


''c/4 - T 


-A 2 
w 

8 cos A 


6 X sin A 


c/4 


c A 
w w 


^ 1 ] 


per deg 


(5. 2. 1.1) 


where, as obtained from Table 2.1: 
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A is the wing aspect ratio 
w 

is the sweep of the wing quarter-chord line 


c is the wing mean aerodynamic chord 
w 

X is the location of the wing aerodynamic center behind the 


center of gravity on the mean aerodynamic chord. 


Since the wing on the S2R-800 is unswept, Equation 5. 2. 1.1 reduces to: 


- 229.2 “ A '*'8 (5.2.1.2) 

6 w 


where: 


C is the wing lift coefficient from Figure 5. 2. 1.1. 

Li 

w 

The contribution of the wing to C of the subject airplane is 

3 

calculated in Table 5. 2. 1.1 to be 5 

(C ) “ .000239 per deg 


5.2.2 Fuselage Contribution to C 

The net contribution of the fuselage and wing-fuselage inter- 
ference to C , based on wing area and wing span and referenced to 
a selected center-of-gravity position, may be obtained from the 
following equation; 


^^f^s *"f 

■ -h — r 


where: 


(S^)g is the fuselage side area, from Figure 2.1.1 

S is the wing area, from Table 2,1 
w 


1 



is the fuselage length, from Figure 2.1.1 
Kjj is an empirical correlating factor for fuselage plus wing- 
fuselage interference, obtained from Figure 5. 2. 2.1. 

The contribution of the fuselage of the subject airplane to C 

n 

I 

is calculated in Table 5. 2. 2.1. 


5.2.3 Vertical Tail Contribution to C 

The contribution of the vertical tail to the weathercock stability 
in the presence of the wing, fuselage, and horizontal tail is obtained 
from: 


^^n„^v(wfh) “^^y.^v(wfh) I 


3 

where; 


a cos a - Z sin a 

V V 

b 

w 


per deg 


(5. 2. 3.1) 


(Cy. the contribution of the vertical tail to the side 

3 

force due to sideslip, obtained from Section 5.1.3 

are the distances from the center of gravity to the quarter 
chord of the vertical tail mean aerodynamic chord, parallel and per- 
pendicular, respectively, to the x-body axis with positive below 
the center of gravity, obtained from Figure 2.1.4. 

The contribution of the vertical tail to the weathercock stability 
of the subject airplane is calculated in Table 5. 2. 3.1. 


5.2.4 Winglet Contribution to C 

Winglet contribution to C is obtained in a similar manner to 

"3 

the vertical tall calculations. 




*WLI 


w 


in a ^ 


per deg (5. 2. 4.1) 
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where : 


(C ) is the side force due to sideslxp of the winglets, 

' WLT 

p 

from Section 5.1.4 

{, z. are the distances from the center of gravity to the 
WLT, 

quarter chord of the winglet mean aerodynamic chord, parallel and 
perpendicular, respectively, to the x-body axis with positive 

below the center of gravity, obtained from Figure 2.2. 

This approximate method does not take into account winglet—wing 

interference or changes in wing span loading that winglets produce. 

Both of these effects can be significant. 

The contribution of the winglets to C of the subject airplane 

6 

is calculated in Table 5. 2. 4.1. 

5.. 2. 5 Weathercock Stability of the Complete Airplane, Winglets 
on and Winglets Off 

The weathercock stability, C , of the subject airplane, winglets 
and power off, is given by Equation 5.2.1. Table 5. 2. 5.1 summarizes 
the calculations and lists the effects of winglet cant angle on . 

V 


5.10 



5.3 Rolling Moment Derivative, C 

e 

The airplane rolling moment due to sideslip, C. , is composed of 
the following contributions: 

(1) Wing 

(2) The effect of the fuselage on the wing contribution 

(3) Vertical tail 

(4) Winglets 

These contributions to C« can be represented as : 

B 

^S.^PROP OFF ” **' (w) ^^i,.^v(wfh) (5.3.1) 

WLT OFF ^ ^ ^ 

^^«,„^PROP OFF “ ^^2, ^ PROP OFF ^^Jl-^WLT (5.3.2) 

^ WLT ON ^ WLT OFF ^ ’ 


5i3.1 Wing Contribution to C. 

6 

At low angles of attack and subsonic speeds, the dihedral effect 
contribution by the wing is primarily a function of wing aspect ratio, 
taper ratio, and dihedral angle. 


‘=1 

(-cT’a * '■« 

B w L w 


(5. 3. 1.1) 


where: 


Cj^ is the wing lift coefficient, from Figure 5. 2. 1.1 


w 


6 

("p — ) is the effect of uniform geometric dihedral on C , obtained 
w 

from Figure 5. 3. 1.1 

% 

(— — ) is the aspect ratio contribution to C„ , obtained from 
^ 3 

Figure 5. 3. 1.2 
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r is the wing geometric dihedral, from Table 2.1« 

The contribution of the wing to C. of the subject airplane is calcu- 

B 

lated in Table 5.3.1.1* 


5,3.2 Effect of Fuselage on Wing Contribution to 


While the contribution of the fuselage alone to C. is negligible, 

8 

the fuselage does influence the flow over the wing which can alter the 
wing contribution significantly. Equation 5. 3. 2.1 from Reference 1 
accounts for this wing-fuselage interference. 




1.2»^ Z , , ^ 

w^ (b + .S ) _ .OOOSv'A ( — r per deg 

57.3 b^ ^ b^ '' 


/ c n o 1 \ . 


where: 

(d ) is the diameter of the equivalent circular fuselage at 
' f w 

the wing, obtained from Figure 2.1.1 

Z is the vertical position of wing below centerline of equivalent 
w 

circular fuselage, from Figure 2.1.1 

h is the height of the fuselage at the wing location, obtained 


from Figure 2.1.1 

j 

w is the width of the fuselage at the wing location, obtained 
from Figure 2.1.1. 

The wing-fuselage interference effects on of the subject airplane 

8 

are calculated in Table 5. 3. 2.1. 

5.3.3 Vertical Tail Contribution to C* 

8 

The vertical tail contributes to the airplane C by virtue of 

8 

the rolling moment produced by the vertical tail side force due to 


5.12 



sideslip. Equation 5. 3.3.1 from Reference 1 Is used to determine 
the vertical tall contribution. 


^v(wfh) 


^Sg^v(wfh)'| 


( Z cos 0 + 1 sin o 

‘ V V 


w 


f 


per deg 


(5. 3. 3.1) 


where; 


(C ) / irt-\ is the vertical tail side force due to sideslip 
v(wfh) 

in the presence of the wing, fuselage, and horizontal tail, from 
Table 5. 1.3.1 

is the perpendicular distance from the x-body axis to the 
quarter chord of the vertical tail mean aerodynamic chord, from 
Figure 2.1.4 

1^ is the distance along the x-body axis from the center of 
gravity to the quarter chord of the vertical tail mean aerodynamic 
chord, from Figure 2.1.4. 


The contribution of the vertical tail to the C. of the subject 
airplane is calculated in Table 5. 3. 3.1. 


5.3.4 Wlnglet Contribution to C. 

B 

The contribution of the winglets to C. is calculated in a 

6 

similar manner as the vertical tall contribution. Since the method 

presented in this section does not account for span loading variations 

induced by the winglets, it can be assumed that the method will tend 

to be inaccurate. Also, separation and interference effects are 

neglected, which can lead to large errors. However, since no analytical 

methods exist for predicting winglet contribution to C. , Equation 

3 

5. 3. 4.1 will be used to approximate it. 
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^WLT ” ^WLT 

P P 


( °°° ° * ‘■ait ° 


V 


w 


s 


per deg (5. 3. 4.1) 


where; 


(C ) is the winglet side force due to sideslxp, obtained from 
'• y. WLT 

p 

Table 5. 1.5.1 

2^^ is the perpendicular distance from the x-body axis to the 

quarter chord of the winglet mean aerodynamic chord, from Figure 2.2 

A is the distance along the x-body axis from the center of 
WLT 

gravity to the quarter chord of the winglet mean aerodynamic chord, 
from Figure 2.2. 


The contribution of the vertical tail to the 
plane is calculated in Table 5. 3. 4.1. 


6 


of the subject air- . 


5.3.5 C of the Complete Airplane, Winglets On and Winglets Off 

The dihedral effect, C. , of the complete airplane, winglets 

B 

and power off, is given by Equation 5.3.1. Table 5. 3. 5.1 summarizes 
the calculations and lists the effects of winglet cant angle on C, 

B 

of the subject airplane. 
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* 

Table 5. 1.2.1: Fuselage Contribution to G 

^3 


Symbol 

Description 

Reference 

Magnitude 

z 

w 

Distance from body centerline to 1/4 point 
of exposed wing root chord, m (ft) 

Figure 2.1.1 

.43 (1.4) 

h 

Maximum body height at wing-body inter- 
section, ra (ft) 

Figure 5.2.2. I 

1.62 (5.32) 

K. 

1 

Wing fuselage interference factor 

Figure 5. 1.2.1 

1.25 

V 

Fuselage volume, (ft^) 

Table 2.1 

7.72 (254.8) 

S 

w 

Wing area, m^ (ft^) 

Table 2.1 

30.34 (326.6) 

(c ) 
y/B 

Body side force due to sideslip, 
per deg 

Potential 

theory 

-0.0195 



- • 

Table 5. 1.3.1: Vertical Tail Contribution to C 

^0 

(a) Effective Aspect Ratio 


Symbol 

Description 

Reference 

Magnitude 

Sh 

Horizontal tail area, (ft^) 

Table 2.1 

5.25 (56.47) 

s 

V 

Vertical tail area, (ft^) 

Table 2.1 

2.12 (27.77) 

b 

V 

Vertical tall span, m (ft) 

Table 2.1. 

1.51 (4.96) 


Vertical tall aspect ratio 

Table 2.1. 

1.08 

tA 

Vertical tall chord at horizontal tail, 
m (ft) 

Figure 2.1.4 

1.78 (5.84) 


Distance from leading edge of vertical tail 
to aeCe of horizontal tail, in plane of 
horizontal tall, m (ft) 

Figure 2.1.4 

.3048 (1.0) 


Distance from root chord of vertical tail 
to root chord of horizontal tail, m (ft) 

Figure 2.14 

-3048 (-1.0) 

Wj)v 

Depth of fuselage at quarter-root chord 
of vertical tail, m (ft) 

Figure 2. 1.1 

.61 (2.0) 

A^(f) 

Ratio of vertical tail aspect ratio In 
presence of fuselage to isolated vertical 
tail aspect ratio 

Figure 5. 1.3.1 
(a) 

1.53 

A^(£h) 

A^(f) 

Ratio of vertical tall aspect ratio in 
presence of fuselage and horizontal tail 
to aspect ratio of tall in presence of 
fuselage alone 

Figure 5. 1.3.1 
(b) 

0.925 

•s« 

Relative tail size factor 

Figure 5. 1.3.1 

1.18 


Summary; A - 1,5I 
^'eff 


5.16 



Table 5. 1.3.1: 


(Continued) 

(b) Vertical Tail Lift Curve Slope 


Symbol 


Description 

Reference 

Magnitude 

A 

Effective vertical tail aspect ratio 

Table 

2.1 

1.51 

^eff 







Vertical tall 

half-chord sweepback, deg 

Table 

2.1 

-3.4 

o 

Vertical. tail 
per rad 

section lift curve slope, 

Table 

2.1 

5.655 

k 

V 

a 

2tt 




0.9 

Summary: 

(C, ) V “ 2.145 

L v(fh) 
a 

per rad “ 0.0374 per deg 






5.18 


Table 5. 1.3.1: (Concluded) 

(c) Vertical Tail Contribution to C 

^8 


Symbol 


Description 

Reference 

Magnitude 



Wing reference-area, m^ (ft^) 

Table 2.1 

30.34 (326.6) 

<^M>» 


Vertical tail quarter-chord sweep angle, 
deg 

Table 2.1 

14 

z 

w 


Distance from equivalent fuselage 
centerline^ m (ft) 

Figure 2.1.1 

0.43 (1.4) 

(Wf>w 


Width of equivalent fuselage at wing, 
m (ft) 

Figure 2.1.1 

1.26 (4.15) 

A 

w 


Wing aspect ratio 

Table 2.1 

5.81 

^ If) 


Wing wake and fuselage sldewash factor 

Equation 5. 1.3. 3 

1.02 

k! 


Relative body size to tail size parameter 

Figure 5. 1.3. A 

0.850 


Sununary: (C ^ -.00225 per deg 
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Table 5. 1.4.1: Wlnglet Contribution to C 

(a) Zero Cant Angle 


Symbol 

Description 

Reference 

Magnitude 

b» 

V 

Tail span above wing plane, m (ft) 

Table 2.1 

1.52 (4.98) 

b 

V 

Total tail span, m (ft) 

Table 2.1 

1.52 (4.98) 

b 

w 

Wing span,, m (ft) 

Table 2.1 

13.27 (43.55) 


Fuselage diameter at wlnglet quarter- 
root chord, m (ft) 

Figure 2,1.1 

1.62 (5.32) 


Length of fuselage, m (ft) 

Figure 2.1.1 

8.5 (27.88) 

^l^LT 

Winglet area, m^ (ft^) 

Table 2.1 

1.8 (19.34) 

A 

V 

Tail (wtnglet) aspect ratio 

Table 2.1. 

1.28 

A 

^eff 

Twin tails (wlnglets) effective 
aspect ratio 

Figure 5. 1.4.1 

1.51 

(c ) 

^B ^eff 

Lift curve slope of one vertical tail 
panel (wlnglet), per rad 

Figure 5. 1.4.1 

2.6 

^Sg^v(wfh) 

(c ) 

^eff 

Mutual interference factor 

Figure 5. 1.4.1 

1.0 


Summary: (C ) = -0.00537 per deg 

CANT=0 
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Table 5. 1.4.1 


Cant Angle , deg 

+20 

0 

-10 

Table 5. 1.4.1 
Cant Angle » deg 


+20 

0 

-10 


+ 




(C , per deg (Equation 5. 1.4. 2) 

y WLT 

^ CANT 


-.000237 

0 

.000118 



(Concluded) 
(c) Summary 



) 

e 


WLT 


^Sg^WLT ■** ^^y AjLT ’ 

^WLT=0 ’^WLT 


per deg 


-0.00561 

-.00537 


-.00525 
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Table 5. 1.5.1: C of the Complete Airplane 

^3 


Symbol 

Description 

Reference 

Magnitude 

(C ) 

Wing contribution, per deg 

Section 5.1.1 

-.00035 

(C ). 

Fuselage contribution, per deg 

Table 

5. 1.2.1 

-.003 

^Sg^v(wfh) 

Contribution of vertical tall In 
presence of wing, body and horizontal 
tall, per deg 

Table 

5. 1.3.1 

-.00225 

^WLT 

^3 ”^^20“ 

Contribution of wlnglets with 20-degree 
cant, per deg 

Table 

5. 1.4.1 

-.00561 

^WLT 

yp wbioo 

Contribution of wlnglets with no cant, 
per deg 

Table 

5. 1.4.1 

-.00537 

^WLT 

^3 -10“ 

Contribution of wlnglets with -10“ cant, 
per deg 

Table 

5. 1.4.1 

-.00525 

Summary: Wlnglets off C = -.0056 per deg 

^3 





Wlnglets on 

20" cant C = -.0112 

^3 

0" C = -.0110 

^3 


- 10 “ 


C = -.0108 
^3 




Table 

5. 2. 1.1: Wing Contribution to C 

8 

Symbol 

Description 

Reference Magnitude 

A 

w 

Wing aspect ratio 


Table 2.1 5.81 

^c/A 

Sweep of wing quarter-chord line, deg 

Table 2.1 0 

c 

Wing mean aerodynamic chord, m (ft) 

Table 2.1 2.29 (7.5) 

Cl 

Wing lift coefficient 

Figure 5. 2. 1.1 f (a) 

w 




Summary: 

<C ) - 0.000239 C 2 

"b " '•w 

per deg 



© 

© 

© 



Figure 5. 2. 1.1 



a, deg 

w 

(C ) = 0.00023902 

“e W 


-A 

-0.025 

.00000598 


0 

0.27 

.00017A 


A 

0.56 

.0000750 


8 

0.86 

.000177 


12 

1.15 

.000316 




Table 5. 2. 2.1: Fuselage Contribution to C 

“3 

Symbol 

Description 

Reference 


Fuselage side area, m^ (ft^) 

Figure 2.1.1 

s 

w 

Wing area, (ft^) 

Table 2.1 


Length of fuselage, m (ft) 

Figure 2.1.1 


Wing span, m (ft) 

Table 2.1 

Z 

w 

Vertical position of wing below center- 
line of equivalent fuselage, m (ft) 

Figure 2.1.1 

(wj 

Width of equivalent fuselage at the 
wing, m (ft) 

Figure 2.1.1 


h 2 Geometric fuselage parameters, m (ft) 

Figure 5. 2. 2.1 

•Si 

Empirical factor for fuselage In 

presence of wing 

Figure 5. 2. 2.1 

Summary : 

(C^ ^f(w) -0.000423 per deg 

3 



Magnitude 

8.36 (90) 
30.34 (326.6) 
8.5 (27.88) 
13.27 (43.55) 
0.43 (1.41) 

1.26 (4.15) 

as listed 
.0024 



5.24 



Table 5. 2, 3.1: Vertical Tail Contribution to C 


Symbol 

Description 

Reference 

Magnitude 

^Sp^vCwfh) 

Contribution of vertical tail to side force 
due to sideslip, per deg 

Table 5. 1.3.1 

-.00225 

A 

V 

Distance along x-body axis from center of 
gravity to quarter chord of vertical tail 
mean aerodynamic chord, m (ft) 

Figure 2.1.4 

5.09 (16.7) 

Z 

V 

Perpendicular distance from x-body axis to 
quarter chord of vertical tail mean aero- 
dynamic chord, m (ft) 

Figure 2.1.4 

-1.07 (-3.5) 

b 

w 

Wing span, m (ft) 

Table 2.1.2 

13.27 (43.55) 

Summary : (C^ 

) , “ 0.000863 cos a + 0.000181 sin a 

p'v(wfh) 

per deg 



® 

a, deg 

(D 

cos^^ 

® 

sln^^ 

^ +0.0001810 

-4 

0.9976 

-0.06976 

.000848 

0 

1 

0 

.000863 

4 

0.9976 

0.06976 

.000874 

8 

0.9903 

0.1392 

.000880 

12 

0.9781 

0.2079 

.000882 
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Table 5.2.A.1: Winglet Contribution to C 

n 


Symbol 

Description 


Reference 

Magnitude 

P 

Contribution of winglets to 
due to sideslip, per deg 

side force 

Table 5. 1.4.1 

As listed 

^WLT 

Distance along x-body axis from center 
of gravity to quarter chord of winglet 
mean aerodynamic chord, ra (ft) 

Figure 2.2 

.44 (1.44) 

^WLT 

Perpendicular distance from 
to quarter chord of winglet 
dynamic chord, m (ft) 

x-body axis 
mean aero- 

Figure 2.2 

-.43 (-1.4) 


Wing span, m (ft) 


Table 2.1.2 

13.27 (43.55) 

Summary : (C^ 

= -(C ),„„ [.0331 cos 

' y„ WLT '■ 

p p 

a + .032 sin a] 

per deg 


© 

© 

0 

© 

© 

© 

o, deg 

cos sin 

^^ri-^WLT 20® 

P 

^^n.^WLT 0“ 

P 

^^n ^WLT -: 

P 

-4 

0.9976 -0.06976 

-.000172 

.000166 

.000162 

0 

1 0 

0.000185 

.000178 

.000174 

4 

0.9976 0.06976 

0.000197 

.000190 

.000185 

8 

0.9903 0.1392 

0.000208 

.000200 

.000196 


0.9781 


12 


0.2079 


0.000218 


000210 


000205 



Table 5. 2. 5.1: C of the Complete Airplane; 


, deg 

^PROP OFF 
^ WLT OFF, 

^PROP OFF 
^ WLT ON 20 

per deg 

per deg 

-4 

0.000419 

0.000591 

0 

0.000457 

0.000642 

4 

0.000526 

0.000723 

8 

0.000634 

0.000842 

12 

0.000775 

0.000993 


Winglets On and Winglets 

Off 

^PROP OFF 
^ WI.T ON 0", 
per deg 

^PROP OFF 
^ WLT -10° , 
per deg 

0.000585 

0.000581 

0.000635 

0.000631 

0.000716 

0.000711 

0.000834 

0.000830 

0.000985 

0.000980 



Table 5. 3. 1.1: Wing Contribution to 


Symbol 

Description 

Reference 

Magnitude 


Wing lift coefficient 

Figure 5. 2. 2.1 

f (a) 


Wing aspect ratio 

Table 2.1. 

5.81 

X 

V 

Wing taper ratio 

Table 2.1 

1.0 

^/2 

Sweep of wing half-chord line, deg 

Table 2.1 

0 


Wing geometric dihedral, deg 

Table 2.1 

3.5 

( 

V 

Low-speed variation of C. as a 

6 

function of C- , per deg 
^w. 

Figure 5. 3. 1.1 

-.0018 

% 

Effect of geometric dihedral on , 

per (deg)^ 

Figure 5. 3. 1.2 

-.00021 

Summary: 

(C ) - -0.000735 - .0018 C. per deg 

^3 ^ w 




© 


© 

a, deg 

V 

(C, )- -0.000735 + 

' *’3 w 

^ -0.0018 ® 

-4 

-0.025 

-0.00069 

0 

0.27 

-0.00122 

4 

0.56 

-0.00174 

8 

0.86 

-0.00229 

12 

1.15 

-0.00280 
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Table 5. 3. 2,1: Effect of Fuselage on Wing 

Contribution to C. 

e 


Symbol 

Description 

Reference 

Magnitude 

A 

w 

Wing aspect ratio 

Table 2.1 

5.81 

b 

w 

Wing span, m (ft) 

Table 2.1 

13.27 (43.55) 

Z 

w 

Vertical position of wing below centerline 
of equivalent circular fuselage, m (ft) 

Figure 2.1.1 

0.43 (1.4) 

(dp) =h=w 
f w 

Diameter of equivalent circular fuselage 
at wing, m (ft) 

Figure 2.1.1 

1.26 (4.15) 

r 

w 

Wing geometric dihedral, deg 

Table 2.1 

3.5 


Summary: 




= 0,000271 per deg 
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Symbol 


Table 5. 3. 3.1: Vertical Tail Contribution to C* 

e 

Description Reference Magnitude 


^^yg^v(wfh) 

Vertical side force due to sideslip in 
presence of wing, fuselage, and horizontal 
tail, per deg 

Table 5. 1.3.1 

-0.00225 

z 

V 

Distance from x-body axis to quarter chord 
of vertical mean aerodynamic chord, ra (ft) 

Figure 2.1.4 

-1.07 

(-3.5) 

z 

V 

Distance along x-body axis from center of 
gravity to quarter chord of vertical-tail 
mean aerodynamic chord, m (ft) 

Figure 2.1.4 

5.09 

(16.7) 

b 

w 

Wing span, m (ft) 

Table 2.1 

13.27 

(43.55) 


Summary : 


^^£.^v(wfh) 

p 


-0.000181 cos a + .000863 sin a, per deg 


© 

a, deg 


© ® 

cos 0 sin 0 (C^ = -0.000181 0 + 

8 (vfh) +0,000863 0 


-4 

0.9976 

-0.06926 

-0.000241 

0 

1 

0 

-0.000181 

4 

0.9976 

0.06976 

-0.000120 

8 

0.9903 

0.1392 

-0.000059 

12 

0.9781 

0.2079 

0.000002 


i % 
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Table 5. 3. 4.1: Winglet Contribution to C. 

6 


Symbol 


Description 


Reference 

Maenitude 

^So^WLT 

P 


Contribution of winglets to 
due to sideslip, per deg 

side force 

Table 5. 1.4.1 

As listed 

*'WLT 


Distance along x-body axis from center 
of gravity to quarter chord of winglet 
mean aerodynamic chord, m (ft) 

Figure 2.2 

.44 (1.44) 

^WLT 


Perpendicular distance from 
to quarter chord of winglet 
dynamic chord, m (ft) 

x-body axis 
mean aero- 

Figure 2.2 

-.43 (-1.4) 



Wing span, m (ft) 


Table 2.1 

13.27 (43.55) 

Summary : 


>WLT “ ^WLT t“0-0321 cos a + 0.0331 sin 

3 3 

a] per deg 


© 

© 

0 

© 

© 

© 


deg 

cos sin 

20“ 

P 

^^«,-^WLT 0“ 

P 

^^il.^WLT -10“ 
3 


-4 

0.9976 -0.06976 

-0.000192 

-0.000184 

-0.000180 


0 

1 0 

-0.000180 

-0.000172 

-0.000168 


4 

0.9976 0.06976 

-0.000167 

-0.000159 

-0.000155 


8 

0.9903 0.1392 

-0.000154 

-0.000146 

-0.000142 


12 

0.9781 0.2079 

-0.000137 

-0.000131 

-0.000128 
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Table 

5. 3.5.1: C. of 

3 

the Complete Airplane: 

# 

Wlnglets On and Wlnglets 

Off 

a, deg 

‘^PROP OFF 
WLT OFF, 
per deg 

'‘a 

'^PROP OFF 
WLT 20® , 
per deg 

’^PROP OFF 
WLT 0®, 
per deg 

"^PROP OFF 
WLT -10® 
per deg 

-4 

-0.00066 

-0.000852 

-0.000844 

-0.000841 

0 

-0.00113 

-0.00131 

-0.00130 

-0.00129 

4 

-0.00160 

-0.00178 

-0.00176 

-0.00176 

8 

-0.00208 

-0.00223 

-0.00222 

-0.00222 

12 

-0.00253 

-0.00267 

-0.00266 

-0.00265 




a 

























For subject airplane - 



Figure 5. 2. 2.1: Empirical factor related to the derivative 

C for the fuselage plus wing-fuselage inter- 
ne 

f erence (Reference 1) . 



prop off 
vlt off 


O O 


DATCOM 


WIND TUNNEL 


Figure 5. 2.5.1! 


a - deg 

Comparison of predicted C » 

winglets off, to full scale wind 
tunnel data for the Ayres Thrush 
ROR_flnn . 







prop off 
wit off 



DATCOM 

per deg 

-.OOU 

-.005 

J_ 

(3 WIND TUNNEL 
1 1 1 


-1 

^ 0 

U 8 12 1 

a - deg 


Figure 5.3. 5.1: 

Comparison of predicted C» , 

6 

winglets off, to full scale wing 
txmnel data for the Ayres Thrush 
S2R-800. 







CHAPTER 6 


PREDICTION OF LIFT AND STABILITY CHARACTERISTICS 
USING A COMPUTER METHOD 

In this chapter the Quasi-Vortex Lattice Method (QVLll) , developed 

by Lan (Reference 3) is used to predict the power-off lift curve slope 

and zero lift angle for the S2R-800. C , C and C. are also calcu- 

B B 

lated. 

6.1 Description of QVLM 

In QVLM the lifting surface is divided into a number of small 
lifting elements. The continuous vortex distribution representing 
the wing in a uniform flow is replaced with a discrete one. Wing 
edge square— root singularities and the Cauchy singularity in the 
downwash integral are theoretically accounted for. A mathematical 
description of QVLM is included in Reference 3. 

6.1.1 Program Capabilities 

The QVLM program used to predict stability derivatives for the 
Ayres Thrush has these following noteworthy features: 

1) It is applicable to nonplanar wing configurations, such 
as wing-wlnglet combinations. 

2) It cannot account for the fuselage effect on the 
stability derivatives. 

3) If the airplane tails do not have camber, their effect 
on stability derivatives cannot be computed. 

A) Arbitrary wing camber shapes defined at three spanwise 
stations or less are used in the program through cubic 


6.1 


spine interpolation. 

5) The vortex-lift effect is calculated through the use of 
Polhamus' suction analogy. 

6.2 QVLM-Predicted and Oq 

a 

For computing , QVLM assumes attached potential flow and 
J-i 

a 

sets the wing angle of attack to 1 radian. The program outputs 

where • Since attached potential flow is assumed, separation 

L li 

a 

effects and stall behavior cannot be predicted. For the Ayres Thrush 
C,. = 0.0715 per deg as predicted by QVL!! (Reference 3). 

For computing Oq, several C^^'s were computed at different 
angles of attack. The zero lift angle could then be found by inter- 
polation. For the Ayres Thrush; 

^“O^QVLIi “ 

These results are compared to full-scale wind-tunnel data in 
Figure 6.2.1. 

6.3 QVLM-Predicted C , C , and C. 

^6 "S 6 

The QVLM calculations were done for winglets on and winglets off. 
The winglet cant angles were +20°, 0®, and -10®. Since the version 
of QVLM used for these calculations could not account for fuselage 
effects, it was not possible to compare these results with the 
available S2R-800 wind-tunnel data. The results are plotted 'in Figures 
6.3.1 to, 6. 3. 3. It is noted that Reference 4 (Figure 39) contains 
small scale model wind tunnel data on the effect of winglet cant on C 
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-h 0 4 8 12 l6 20 


a - deg 

Figtire 6.2.1: Lift curve of the S2R-800 based on 

QVLM computation. 
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J 


WING ALONE 


.20 deg cant 
0 deg cant 
-10 deg cant 

0 deg 



a - deg 

Figiire 6.3.1: Effect of winglet cant angle on C based 

yg 

on QVLM computation for the S2R-800. 


WING + WLT 

• — WING + WLT 
WING + WLT 

I ^ 
'WLT 
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WING ALONE 


WING + WLT 20 deg cant 
WING + WLT 0 deg cant 
WING + WLT -10 deg cant 

'wLT “ ® ‘‘'3 



-2 0 2 k 6 

a - deg 

Figure 6.3.2: Effect of winglet cant angle on 


on QVLM computation for the S2R-800. 


WING ALONE 


20 deg cant 
0 deg cant 


•t 


per deg 


WING + WLT 

— “ — WING + WLT 

WING + WLT -10 deg cant 

'hut ■ ” ''=3 



a - deg 

Figure 6 . 3 . 3 J Effect of winglet cant angle on C 

based on QVLM computation for the 
S2R-800. 
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CHAPTER 7 


CONCLUSIONS AND RECOMMENDATIONS 


In this report an analytical method (Reference 1) and a computer 
method (Reference 3) were presented and used to calculate selected 
lift characteristics and sideslip derivatives of the Ayres Thrush 
S2R-800. Where possible, these results were compared to full-scale 
wind-tunnel data and to each other. Based on these comparisons, the 
following conclusions and recommendations are made: 

1. The analytical method of Reference 1 and computer method 

of Reference 3 both give C^ and Oq values which correlated well 

a 

with the wind-tunnel data. The analytical method gives a good "ball- 
park" estimation of the nonlinear lift behavior, while the computer 
method of Reference 3 does not apply to the nonlinear region. 

2. The analytical predictions of Reference 1 for the sideslip 

derivatives C , C and C, show good agreement with wind-tunnel 
^6 6 B 

data. The predicted C is within 10 percent of the average tunnel 

C value. While the trends with angle of attack are opposite, the 
^B 


nominal predicted agrees with the average tunnel . Separation 

B B 

and interference effects caused by the fuselage need to be incorporated 

in the analytical method. The predicted C. compares well with wind- 

B 

tunnel results. 

3. The predictions by the method of Reference 3 indicate that 
winglet cant angle has a significant influence on the sideslip deriva- 
tives. From Figures 6. 3. 1-6. 3. 3 it is evident that the inwardly canted 
winglets have the least effect on the lateral-directional stability of the 
Ayres Thrush. Therefore, it appears that the use of inwardly • 
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canted winglets may offer a means of minimizing chemical spray wake 
interaction effects without degrading handling qualities (see 
Reference 4). 

4. Computer results show that the analytical method of ref- 
erence 1 is deficient in predicting winglet effect on the sideslip 
derivatives. It is recommended that a more accurate analytical method 
for calculating winglet effect on the sideslip derivatives, incor- 
porating the following procedure, be developed. 

a) To account for wing span loading increase due to 
winglets (endplate effect), calculate a modified wing aspect 
ratio based on winglet area. 

b) Calculate winglet C^^^ based on airfoil section 

a 

properties and adjusted for finite span. 

c) Based on wind-tunnel data, calculate and plot wing- 
winglet interference factors. 

d) Using the values calculated above, follow the 
procedure developed in this report, making modifications 
where common sense dictates. 
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